Entangling many atomic ensembles through laser manipulation 
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We propose an experimentally feasible scheme to gener- 
ate Greenberger-Horne-Zeilinger (GHZ) type of maximal en- 
tanglement between many atomic ensembles based on laser 
manipulation and single-photon detection. The scheme, with 
inherent fault tolerance to the dominant noise and efficient 
scaling of the efficiency with the number of ensembles, al- 
lows to maximally entangle many atomic ensemble within the 
reach of current technology. Such a maximum entanglement 
of many ensembles has wide applications in demonstration of 
quantum nonlocality, high-precision spectroscopy, and quan- 
tum information processing. 

PACS numbers: 03.65.Ud, 03.67.-a, 42.50.Gy 42.50.-p 

Quantum entanglement links two or more distant sub- 
systems in a profound quantum mechanical way. Such 
a link has found wide applications in demonstration of 
quantum nonlocality ||lj,2|, high-precision spectroscopy 

, and quantum information processing including com- 
putation, communication and cryptography There 
are great experimental efforts recently to get more and 
more subsystems entangled j^-^, since with more sub- 
systems entangled, quantum nonlocality becomes more 
striking QJ^, and the entanglement is more useful for 
various applications [^-^. In most of the experimental 
efforts, the subsystems are take as single-particles, and 
up to now three to four atoms or photons have been en- 
tangled with a linear ion-trap ||^, with a spontaneous 
parametric down converter |^,^, or with a high-Q cav- 
ity §. There are also proposals to entangle indistin- 
guishable atoms in Bose- Einstein condensates , or to 
weakly entangle two macroscopic atomic ensembles [0, 
and the latter has been demonstrated in a recent exciting 
experiment |12|. 

In all the experimental efforts, it is hard to continu- 
ously increase the number of the entangled subsystems 
due to the fast exponential decrease of the preparation 
efficiency |6|j^,^ or due to noise and imperfections in the 
setup 0,^. Here, we propose a scheme to generate GHZ 
type of maximal entanglement between many atomic en- 
sembles with the following features: firstly, the scheme 
has built-in fault-tolerance and is robust to realistic noise 
and imperfections. As a result, the physical requirements 
of the scheme are moderate and well fit the experimen- 
tal technique. Secondly, the preparation efficiency of the 
GHZ entanglement only decreases with the number of en- 
sembles by a slow polynomial law. Such an efficient scal- 



ing makes it possible to maximally entangle many (such 
as tens of) ensembles with the current technology. Our 
scheme is based on Raman type laser manipulation of the 
ensembles and single-photon detection which postselects 
the desired entangled state in a probabilistic fashion. In 
contrast to the belief that entangling schemes based on 
postselections will necessarily suffer from the fast expo- 
nential degradation of the efficiency, we design a scheme 
which circumvents this problem by making use of quan- 
tum memory available in atomic internal levels. 

The basic element of our system is an ensemble of 
many identical alkali atoms, whose experimental realiza- 
tion can be either a room-temperature atomic gas 13 
or a sample of cold trapped atoms 0,^. The rele- 
vant level structure of the atom is shown by Fig. 1. 
From the three levels \g) ,\h) ,\v), we can define two col- 
lective atomic operators s = (l/\^Na) X^i^i \9)i i^l '^ith 
s = h,v, where 3> 1 is the total atom number. The 
atoms are initially prepared through optical pumping to 
the ground state \g), which is effectively a vacuum state 
|vac) of the operators h, v. The h, v behave like indepen- 
dent bosonic mode operators as long as most of the atoms 
remain in the state \g) . A basis of the "polarization" 
qubit (in analogy to the language for photons) can be de- 
fined from the states \H) = |vac) and \V) = v'^ |vac), 
which have an experimentally demonstrated long coher- 
ence time [p^-p^. Single-bit rotations in this basis can 
be done with high precision by shining Raman pulses or 
radio-frequency pulses on all the atoms. The excitations 
in the mode h can be transferred to optical excitations 
Jl6t and then detected by single-photon detectors. Such a 
transfer has a high efficiency even for a free-space ensem- 
ble due to the collectively enhanced coherent interaction 
as has been demonstrated both in theory ||l6| and in 
experiments 13,1^]. 

The first step for generation of many-party entan- 
glement is to share an excitation between the modes 
hi,hj in two distant ensembles This can be read- 
ily done through a scheme in the recent quantum re- 
peater proposal [T^ , where one prepares the state 

(hi + e"^^^h]^ /V2 |vac) [@, with (j)^^ = <j)j - cj),, an un- 
known phase difference fixed by the optical channel con- 
necting the i,j ensembles. This state, after a single-bit 
rotation, can be transferred to the useful form 
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The basic idea of the preparation scheme in Ref. ||T^ is 
as follows: one excites the ensembles i and j respectively 
through a short weak Raman pulse applied to the transi- 



1 



tion \g) — > |e) so that the forward scattered Stokes hght 
from the transition |e) \h) has a mean photon num- 
ber much smaller than 1. The forward-scattered Stokes 
lights from the two ensembles are then interfered at a 
beam splitter and further detected by two single pho- 
ton detectors. If we successfully get a detector click, we 
do not know from which ensemble this registered pho- 
ton comes from, and due to this indistinguishability, the 
accompany collective atomic excitation should be dis- 
tributed over the ensembles i and j with an equal prob- 
ability amplitude, and we get the state described before. 
This preparation scheme has the following two features: 
first, the preparation only succeeds with a controllable 
small probability po for each Raman driving pulse, and 
needs to be repeated in average l/po times for the final 
successful state generation, with the total preparation 
time to ^ 1/ (pofp), where fp is the repetition frequency 
of the Raman pulses. Second, the scheme, with inher- 
ent resilience to noise, is well based on the current tech- 
nology of laser manipulation. We can safely use it as 
our first step, to generate the state (1) with a fidelity 
F = 1 — Po very close to the unity by controlling the 
probability pq. For instance, with a typical repetition 
frequency fp = lOMHz, one may prepare the state (1) 
with a fidelity F = 1 — po ~ 99% in a time to ~ lO^s. 

Based on the preparation of the state (1), now we show 
how to generate effective many-party entanglement be- 
tween n such atomic ensembles. We prepare the state 
(1) between the i and i -I- 1 ensembles for each i from 1 
to n, and get the following state 

n 

Ivf) = (l/^/2^) n (/it + e^<*-+^4+i) |vac) , (2) 

i=l 

where for convenience we have assumed the notation 
n -|- 1 = 1 for the subscripts, and have used the same 
symbol |vac) to denote the vacuum of the whole n en- 
sembles. In the expansion of the state (1), there are only 
two components which have one excitation on each en- 
semble. This component state is given by 

(n n \ 

l[''l + e^^^l[vl\\ve.c) (3) 
1=1 i=l / 

with (j)t = X]"=i (f'i.i+ij which is exactly the n-party GHZ 
type maximally entangled state in the polarization basis. 
Note that for any practical application of the GHZ entan- 
glement [Q-Q , the state preparation should be succeeded 
by a measurement of the polarization of the excitation 
on each ensemble, which can be done for our system by 
combining single-bit rotations, such as Hadamard trans- 
formations, with the number detection of the mode hi 
through single-photon detectors. If in this measurement 
we only keep the results for which an excitation appears 
on each ensemble (i.e., postselect the case when the de- 
tector on each side registers a click), the states (2) and 
(3) become effectively equivalent since the other compo- 
nents in the state (2) have no contributions to the mea- 
surement. Through this postselection technique, we can 



simply prepare the state (2), which, whenever we put 
it into applications, yields effectively the GHZ entangle- 
ment described by the effective state (3). Here and in the 
following, we call a component of the full state as the ef- 
fective state if only this component has contributions to 
the application measurements (which are the measure- 
ments required for the detection or application of the 
generated state). The effective state is the state postse- 
lected by the application measurements. 

For applications of the GHZ entanglement, we need 
also to know the phase 4>t in the effective state (3), which 
is fixed by the whole setup and in principle can be mea- 
sured. However, a better way is to directly cancel this un- 
known phase 4>t with the following method. Assume that 
we have an even number n of the ensembles. The pair of 
ensembles i and i' = n + 2 — i are put in near proximity so 
that the ensembles i, i+1 and i', {i + 1)' can be connected 
through the same optical channel, which fixes the phases 
to satisfy the relation = <l>i',{i+iy = -</'(j+i)',i' 

With this relation, the accumulated phase (j)t is exactly 
canceled to zero. 

The above preparation scheme of the effective GHZ en- 
tanglement is robust to realistic noise and imperfections. 
The dominant noise in this system is the photon detec- 
tor inefficiency, the transferring inefficiency (induced by 
the spontaneous emission loss) of the excitation from the 
atomic mode hi to the optical mode, and the small de- 
cay of the atomic excitation in each ensemble. All the 
above noise is well described by loss of excitations with 
a overall loss probability denoted by 77. Note that by in- 
cluding the detector inefficiency, we have automatically 
taken into account that the single-photon detectors can- 
not perfectly distinguish between single and two photons. 
It is easy to see that loss of excitations only has influence 
on the success probability to register an excitation from 
each ensemble. Whenever the excitation is registered, its 
polarization is still perfectly entangled as shown by the 
effective state (3). 

Now we consider the efficiency of this scheme, which 
can be described by the total time needed to successfully 
register the effective GHZ entanglement. The prepa- 
ration of the factor state (1) is probabilistic, however, 
due to the available quantum memory provided by the 
metastable atomic modes /i, v, the preparation time ti of 
the state (2) is at most nto if its factor states are prepared 
one after the other, and can be reduced to ti ~ to (in the 
case of n < 1/po) if its factor states are prepared indepen- 
dently at the same time. In contrast to this, in the case of 
no quantum memory, one would need about 1 /po repeats 
of the Raman pulses for a successful preparation of the 
state (2), and a total time to/po~^ ^ to- After prepa- 
ration of the state (2), the projection efficiency (success 
probability) from the state (2) to the effective GHZ state 
(3) is given by (1 — 77)" /2^~^, where we have assumed 
the same loss probability 77 for each ensemble. So the 
total time for registering the 7^-party GHZ entanglement 
is T ~ to2"~^/ (1 — 77)", which increases with the num- 
ber of ensembles exponentially by the factor 2/(1 — 77). 
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Note that this increase has been much slower than the 
case for spontaneous parametric down conversion where 
the exponential increasing factor is about 2 orders larger 
due to the absence of quantum memory j^J^. 

We can in fact further improve the scheme to get a 
much more efRcient scaling of the efficiency, with the time 
T increasing with the party number n only polynomially. 
The improved scheme is divided into the following three 
steps: 

(i) We start with two pairs of ensembles 1,2 and 3,4, 
prepared in the state 1^1,2) ® 1*3,4) with in the 

form of Eq. (1). We then connect these two disjoint pairs 
by preparing the state |^'2,3)- The ensembles 2 and 3 will 
not be involved any more in the following steps for state 
preparation, so we can immediately put them into appli- 
cations by doing the same type of measurements on them 
as if we had generated n-party GHZ entanglement. In 
these measurements, if one excitation is registered from 
each ensemble 2 and 3, we succeed and will go on with 
the next step. Otherwise, we simply repeat the above 
process until we succeed. Upon success, only the com- 
ponent |^'i_4) of the state 1*1,2) 1*3,4) ® 1*2,3) has 
contributions to the measurement with 



I* 



1-4, 



= 1 



7V2) (h\h\ 



2^3 



(4) 



where for simplicity we have neglected the phase (/>i,i+i 
since they will finally cancel each other with the method 
described before. If loss of excitations with a loss proba- 
bility 77 is taken into account for detections on the ensem- 
bles 2,3, a registered click might result from two excita- 
tions, and in this case there will be no excitation in the 
ensembles 1 and 4. So with the loss, upon success of step 
(i), instead of 1*1-4) the effective state of the ensembles 
1-4 is actually described by 



Pl_4 = (1*1-4) (*1-4| + ClPvac) / (1 + Cl) 



(5) 



with the vacuum coefficient ci = 2?7, where pvac stands 
for the vacuum component with no excitation in the un- 
detected ensembles 1 and 4. The probability of a success- 
ful detection on both of the ensembles 2 and 3 is given 
by pi = (1 — rff' (1 -|- 277) /4, which means that we need 
to repeat the process in average 1/pi times for the final 
success of step (i). 

(ii) In step (ii) we further extend the number of en- 
tangled ensembles in the effective state (5). Assume that 
we have applied the method of step (i) in parallel to the 
two disjoint sets of ensembles 1-4 and 5-8, with their ef- 
fective states (each in the form of Eq. (5)) denoted by 
/9i_4 and p5_8, respectively. We connect these two sets by 
first preparing the state 1*4,5) (in the form of Eq. (1)) 
and then putting the ensembles 4, 5 into application mea- 
surements as described in step (i). Upon success of these 
measurements with one excitation registered from each 
ensemble, the postselected state of the ensembles 1-8 is 
effectively described by pi_8 which is similar to Eq. (5), 
but with an increased vacuum coefficient and with 1*1-4) 



replaced by |*i.fe) = (l/^/2) (nti' ^^1 + Otz 4) |vac) , 
(/c = 8). Whenever the measurement fails, we repeat 
the whole state preparation from step (i). The above 
connection process can be continued with the number 
n of effectively entangled ensembles doubled for each 
time of connection. After i times connection, we have 
n = 2*+^. The success probability and the new vac- 
uum coefficient of the ith connection are denoted respec- 
tively by Pi and Ci, which satisfy the following recur- 
sion relations with the previous vacuum coefficient Ci_i 

through - (1 - Jif (1 + 277 -h 2c,_i) / [4 (1 -h c,-if 

and Ci = 2ci_i -|- 277. From these recursion relations, we 
have Ci = 27; (2* — l), which, after substituted into pi, 
yields an explicit expression for the repetition number 
1/pi of the 7th connection. 

(iii) After a desired number n — 2*+^ of the ensem- 
bles have been entangled in the effective state pi-m we 
close the loop in the last step by first preparing the state 
|*n,i) (in the form of Eq. (1)) and then putting the 
last two ensembles 77, 1 into application measurements. 
As usual, we keep the results only when one excitation 
appears from each detected ensemble, and this automat- 
ically eliminates contributions from the vacuum compo- 
nent in the state pi_„ . So the effective state of the whole 
set of ensembles postselected by all the application mea- 
surements is still described by the exact GHZ state (3), 
and the application measurement results should reveal 
perfect GHZ entanglement between the 71 ensembles in 
the polarization degree of freedom. The application mea- 
surements on the ensembles 71, 1 in the last step succeeds 
with a probability = (1 — ilf' / [2 (1 + q)], so the whole 
process needs to be repeated in average \/pi times. 

Now we calculate in this improved scheme how much 
time is needed in total for a successful detection of 
the 77-party GHZ entanglement. This time is given 

by Timp = to/ PiPi]\]^2P3 ' with to, the prepara- 
tion time of the state (1). We consider two limiting 
cases. In the first case with a negligible loss proba- 
bility 77 for each ensemble, we have pi = 2pj = 1/2 
and Timp = 2^'+^io = n'^to/'^, which increases with the 
number 77 of entangled ensembles by the slow quadratic 
law. In the second case with a considerably large loss 
probability 77, the total time Timp is approximated by 



imp 



tn 



27,77/(1-7;)' (77/2) 



which 



increases with 77 faster, but still polynomially (or, more 
accurately, sub-exponentially). The basic reason for the 
improvement from the exponential scaling to the much 
slower polynomial scaling is due to that we have divided 
the whole preparation process into many small steps, 
checking in each step whether the preparation is success- 
ful, and repeating this small step instead of the whole 
process if it fails. 

Finally, we briefly discuss the practical implication of 
this proposal. With the improved scheme, for exam- 
ple, we can generate high-fidelity GHZ entanglement over 
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n = 16 ensembles in a time Timp ~ 50ms with a notable 
loss 77 w 1/3 and a typical choice to ~ 10/xs. With such 
a short preparation time Timp, the noise that we have 
not included, such as the non-stationary phase drift in- 
duced by the pumping lase or by the optical channel, 
is negligible. As long as the number n of the ensem- 
bles is not huge, we can also safely neglect the single-bit 
rotation error (below 10"** with the use of accurate po- 
larization techniques for Zeeman sublevels ||2^ and the 
dark count probability (about 10~^ in a typical detec- 
tion time window 0.1/is) of single-photon detectors. Due 
to the efficient scaling of this scheme, one can use it to 
steadily increase the number of entangled ensembles, and 
it seems reasonable to generate GHZ entanglement over 
tens of ensembles with the current technology. Such an 
extraordinary possibility opens up prospects for many 
exciting experiments and applications. 
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FIG. 1. The relevant atomic level structure with \g), the 
ground state, |e) , the excited state, and \h) , \v) the two 
metastable states (e.g., Zeeman or hyperfine sublevels) for 
storing a qubit of information. The three levels \g) , \e) , 
and \h) can be coupled through a Raman process which is 
useful for measurement of the collective atomic excitation in 
the state \h) [16] and for generating preliminary entanglement 
between two ensembles [17]. 
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